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Calculation of weak matrix elements in domain-wall QCD with the 
DBW2 gauge action 

J. Noaki**" for the RBC Collaboration* 

'^RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 

We report the details of our ongoing quenched calculations of weak matrix elements using the combination 
of domain-wall fermions and the DBW2 gauge action on lattices with a~^ ~ 3 GeV. A strategy to avoid the 
problem of fixed topological charge is introduced in generating gauge configurations. After studying the basic run 
parameters and elemental quantities, we present a preliminary result for the kaon B-parameter {Bk)- 


1. Introduction 

The combination of domain-wall fermions and 
renormalization group (RG) improved gauge ac¬ 
tions has been proven to yield fermions with good 
chiral behavior, allowing progress in the calcula¬ 
tion of weak matrix elements (WME) on the lat¬ 
tice 1^. On the other hand, two examples of 
open issues are the scaling behavior of Bk and 
the effect of the charm quark which contributes 
to K TTTT matrix elements through the quark 
loop contraction ||]. Both of these require the 
calculation of WME, using fermions with good 
chiral behavior, on finer lattices than have been 
employed so far. 

In this article, we present the current results 
from an ongoing quenched numerical simulation 
by the RBC Collaboration. It uses domain- 
wall fermions (DWE) and the DBW2 gauge ac¬ 
tion at a scale a~^ w 3 GeV with physical 
spatial size w 1.6 fm. Though the use of the 
RG-improved DBW2 gauge action improves chi¬ 
ral symmetry, it has been recognized that the 
topological charge of DBW2 gauge configurations 
generated by standard Monte Carlo techniques 
evolves very slowly [Q. This fact motivated 
us to implement a strategy to generate gauge 
configurations whose topological charge is well- 
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Figure 1. The gauge configurations used for spec¬ 
trum (open circle) and Bk (filled circle). Shad¬ 
owed bars indicate data employed in Figure |[ 


distributed |^. 

In the following, we discuss the details of our 
numerical simulation and present preliminary re¬ 
sults for the meson spectrum, residual chiral sym¬ 
metry breaking mass Wres and Bk- 

2. Gauge configuration generation 

On a 24^ x 48 lattice, we generate gauge config¬ 
urations using the Wilson action with /3 = 6.25 in 
steps of 10,000 sweeps and use each of these as the 
initial configuration for a DBW2 evolution with 
P = 1.22 in steps of 1000 sweeps. The status 
of our numerical simulation for 15 independent 
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Table 1 

Simulation parameters for observables. 



spectrum 

Bk 

size 

24^ X 48 

Ms 

1.7 

1.65 

Lg 

8 

10 

mja 

0.02, 0.03, 0.04 

0.008 - 0.040 
in step of 0.008 

# sweep 

every 1000 

every 5000 

#config. 

120 

30 


DBW2 series is illustrated in Figure 0. 

The initial Wilson configurations have « 3 
GeV, the same lattice spacing desired for the 
DBW2 ensemble. The part of the topological 
charge produced by physical, long distance prop¬ 
erties of the Wilson lattices should be very close 
to that of a fully sampled DBW2 ensemble, since 
at these weak couplings such a physical quantity 
is independent of the UV details of the action. 
We rely on this physical argument to allow us to 
use the topological distribution generated by the 
Wilson action as a starting point for our DBW2 
ensemble. With the 15 Wilson configurations dis¬ 
cussed here, we find (Qtop) = —0.60± 3.11. (For 
50 configurations, we have (Qtop) = —0.32 ±3.36, 
with —10 < Qtop ± ±7-) For a few of the DBW2 
evolutions, Qtop changed by one unit, while for 
the longest 30,000 sweep evolution, Qtop did not 
change. Thus, for the DBW2 action, we conclude 
that over-relaxed and heat-bath steps would not 
generate a reasonable sampling of topologies at 
these weak couplings. 

3. Hadron spectrum 

Simulation parameters for the spectrum calcu¬ 
lation are summarized in the middle column of 


Table 2 

Spectrum results for M 5 = 1.7, Lg = 8 . 


mpsimv 

0.612(11), 0.698(9), 0.756(7) 

a-i 

2.89(12) GeV 

TOsa /2 

0.0157(11) 

TTT-res 

1.011(35) MeV (Ms = 1.7, Lg = 8 ) 


0.276(10) MeV (Ms = 1.65, L^ = 10) 



Figure 2. mpga versus evolution in the DBW2 
direction (filled symbols) obtained from averag¬ 
ing over 15 configurations in the vertical direction 
shown in Figure Open symbols represent the 
data from the initial Wilson configurations. 

Table |^. To investigate how our strategy works 
and determine optimal parameters for the calcu¬ 
lation of WME, we studied meson masses and the 
residual quark mass Wres which is induced by the 
explicit chiral symmetry breaking due to finite 
Ls- 

It is important to estimate the equilibration 
time for the DBW2 evolutions, starting from the 
initial Wilson configurations. Figure shows the 
results for the pseudoscalar meson mass mpsa for 
each TO/a and averaged over the configurations in¬ 
dicated in Figure |^. In this figure, one observes 
that about 3000 sweeps are necessary for thermal- 
ization, which is consistent with a similar analy¬ 
sis for TOres- Oui' hadron spectrum determination 
only uses data from 3000 sweeps on. 

Table contains the results from our spectrum 
calculation with M 5 = 1.7 and Lg = 8. The lat¬ 
tice scale a~^ was determined from the rho meson 
mass at the chiral limit. We obtained the strange 
quark mass from the relation TOpg = 2BQmf. 
Further study of TOres as a function of domain- 
wall height M 5 , gave 1.65 as the optimal value. 
Using this, we found TOres < 0.3 MeV for the 
DBW2 action with Lg — 10. This value for Lg 
is used in the calculation of WME. 
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Figure 3. time dependence of Bk for rufa = 
0.016 (upper) and 0.024 (lower). The solid lines 
indicate constant fits. 



Figure 4. Lattice value of Bk as a function of 
m/a. The interpolated point (filled symbol) indi¬ 
cates the value at rufa = msaj2. The solid line 
denotes a fit that is described in the text. 


4. Kaon B parameter 
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was calculated with the parameters summarized 
in the right column of Table |^. Quark prop¬ 
agators from Coulomb gauge fixed spatial wall 
sources located at t = 7 and 41 were used in 
our calculation and each propagator is an average 
of ones with periodic and anti-periodic boundary 
conditions in the time direction. 

One finds reasonable plateaus in t for Bk as 
shown in Figure Figure ^ shows the un¬ 
renormalized value of Bk as a function of m/a 
and a fit to the function BK{mfa) = ^o[l + 
C'm/aln(m/a)] -I- where the value of C 

is taken from chiral perturbation theory 
The data are well fit to this function. An un¬ 
correlated fit has x^/dof = 0.14. The interpo¬ 
lated value Bk = 0.539(21) at the physical point 
m/a = vusaj^ is indicated by a filled square. 
We have not yet calculated the non-perturbative 
renormalization factor for Bk- Employing the 
perturbative one 0, the value of Bk in the MS 
NDR scheme can be estimated. 

Previous works calculated Bk{h = 2 GeV) 
with DWF in similar physical spatial volumes [|[ 


^. While our preliminary result is consistent with 
both of these within roughly one and two stan¬ 
dard deviations, respectively, it is smaller than 
results from larger physical volume ||]. 

We thank RIKEN, BNL and the U.S. DOE for 
providing the facilities essential for the comple¬ 
tion of this work. 
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